Holstein polarons in a strong electric field: delocalized and stretched states 
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The coherent dynamics of a Holstein polaron in strong electric fields is considered under different 
regimes. Using analytical and numerical analysis, we show that even for small hopping constant 
and weak electron-phonon interaction, the original discrete Wannier-Stark (WS) ladder electronic 
states are each replaced by a semi-continuous band if a resonance condition is satisfied between 
the phonon frequency and the ladder spacing. In this regime, the original localized WS states can 
become delocalized, yielding both 'tunneling' and 'stretched' polarons. The transport properties of 
such a system would exhibit a modulation of the phonon replicas in typical tunneling experiments. 
The modulation will reflect the complex spectra with nearly-fractal structure of the semi-continuous 
band. In the off-resonance regime, the WS ladder is strongly deformed, although the states are still 
localized to a degree which depends on the detuning: Both the spacing between the levels in the 
deformed ladder and the localization length of the resulting eigenfunctions can be adjusted by the 
applied electric field. We also discuss the regime beyond small hopping constant and weak coupling, 
and find an interesting mapping to that limit via the Lang-Firsov transformation, which allows one 
to extend the region of validity of the analysis. 
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I. INTRODUCTION 

Quantum electronic transport properties in semicon- 
ductor superlattices in high electric fields have been the 
subject of much interest in recent years. Experimen- 
tal studies have reported interesting phenomena, such 
as Bloch oscillations and Wannier-Stark (WS) ladders.u 
Moreover, a number of theoretical predictions havej-been 
made, including negative differential conductivity^ dy- 
namical localizatianj3 and fractional WS ladders under 
DC and AC fields^ to name a few. 

One important issue for the dynamical behavior of a 
real system is the effect of electron-phonon interactions, 
and extensive research in this field has been reported. 
For example, Ghosh et al. and Dekorsy et al. studied 
coupled Bloch-phonon oscillations U The phonon- assisted 
hopping of an electron on a WS ladder was studied in 
A brilliant variational treatment of inelastic quan- 
tum transport was given in |Q, and anomalies in trans- 
port properties under a resonance condition were studied 
in Q . Govorov et al. studied the optical absorption asso- 
ciated with the resonance of a WSJadder and the optical 
phonon frequency in a system £rE2l A similar and fasci- 
nating system of electron-phonon resonance in magnetic 
fields has also been studied extensivelyO 

Much of the work in this area considers incoherent 
scattering of electrons by phonons, which is the relevant 
picture in systems at high temperature. In this paper, 
however, we report a study of the effects of coherent cou- 
pling of an electron to the phonons of the system, likely 
to be the relevant description at low temperatures. We 
concentrate on how the phonon coupling affects electron 
transport in high electric field, in a situation typically 
achieved in semiconductor superlattices, for instance, but 
also important for electrons in polymer chains,t3 and 



even perhaps reachable in a stack of self-assembled quan- 
tum dots.E3 

Based on a Holstein model,0 typically used to describe 
the small polaron in molecular systems, we analyze the 
spectrum of the system and study its transport proper- 
ties. We use here a non-perturbative description, which 
allows one to elucidate the effects of resonant and non- 
resonant phonon fields on the otherwise localized elec- 
trons residing in a WS ladder, for both weak and strong 
electron-phonon coupling. We find that for small hop- 
ping and coupling constant, an interesting regime results 
when two characteristic energy scales in the system co- 
incide: the one corresponding to the energy spacing be- 
tween the WS levels, and the other corresponding to the 
frequency of the phonons. In this resonant case, the prob- 
lem is one of strong mixing between degenerate states 
with rather different properties. The result is a com- 
plex semi-continuous band structure replacing each of 
the original WS ladder 'rungs', where some of the elec- 
tron states become delocalized, despite the strong elec- 
tric field present. In this case, the phonons interacting 
with the electron provide delocalization, unlike the usual 
situation. Moreover, and in contrast to these extended 
states, 'stretched polarons' are highly degenerate and ex- 
hibit strongly localized electronic components, while the 
phonon component is in fact extended throughout the 
structure. When the system is away from the resonance 
condition, a deformed WS ladder will appear, with very 
interesting substructure. The electronic wave functions 
in this regime are all localized, but with a localization 
length that depends linearly on the degree of detuning. 

In all phonon frequency regimes, the rich dynamical 
behavior of the system is due to (or reflected in) the near 
fractal structure of the spectrum. This, as we will discuss, 
has a direct connection to the Cayley tree structure of 
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the relevant Hilbert space. This structure is of course 
contained in the Hamiltonian and reflected also in the 
structure of the eigenvectors. 

We also extend our studies beyond the weak coupling 
and small hopping constant regime, via a Lang-Firsov 
transformation. This results in an interesting mapping of 
the problem from the strong coupling to the weak limit, 
and allows one to extend the results for a wider parameter 
range. A first report of the resonani_case in the weak 
coupling regime has been published. t3 Here, we give a 
detailed theory of such states, including both resonant 
and non-resonant cases, and present the strong coupling 
interaction regime. 

In what remains of the paper, we give in section II a 
description of the model used and the relevant Hilbert 
space in the problem. Section III contains a complete 
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In the case of strong electric field (or small hopping con- 
stant, 2t/A <C 1) in which we concentrate in this section, 
the Hamiltonian can be simplified as 
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analysis of the resulting spectrum for small hopping con- 
stant and weak phonon coupling. The extension beyond 
that regime is given in section 
presented in section |v|. 



II. DESCRIPTION OF THE MODEL 

We consider the Holstein model,Q which describes an 
electron in a one-dimensional tight-binding lattice, inter- 
acting locally with dispersionless optical phonons. More- 
over, the system is subjected to a strong static electric 
field. The Hamiltonian is then given by 

3 3 

+ w^a+a J +7^c+c J -(a+ + a i ), (1) 

3 3 

where t is the electron hopping constant, lu the phonon 
frequency (or energy, with h — 1), 7 the electron-phonon 
coupling constant, £j = —edEj = — jA the site energy, 
E the electric field, and d is the lattice constant. It is 
known that a series of localized WS states will form in a 
strong field (A = edE ^> t), in the absence of interactions 
between electron and phonons. Within tight-binding the- 
ory, the eigenvalues and eigenfunctions of these WS states 
are, respectively Sj = —edEj, and 



fvl A final discussion is where the effective coupling constant is now 
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where J; is the l-th order Bessel function. The WS states 
are localized states with characteristic length 2t/eE. 

It is helpful to introduce creation and annihilation op- 
erators for WS states as 



d 3 = Ji-j(2t/edE)ci 



(3) 



It is easy to show that this transformation is canonical, 
{df,dj} = 5ij, and that the Hamiltonian can be written 
in terms of dj , d+ as 



A = jt/A. 



(6) 



H = -J2j^d+d ]+ J2 



From (§), one can see phonon-assisted hopping between 
the WS states quite clearly, so that in fact the phonons 
introduce derealization of the WS electrons. The Hamil- 
tonian (^) is similar to that in || , but there are some dif- 
ferences. In the current model, when an electron jumps 
from one site to the next, it can not only emit (or ab- 
sorb) a phonon on (from) its site, but also onto (from) 
the nearby site. This more natural description also pro- 
duces a large difference on the relevant Hilbert spaces. 
The dimension of the relevant Hilbert space in our case 
is 2" — 1, where n is the number of local sites in the chain, 
while the size is n for the model in J9|. As we will see, the 
resulting energy spectrum and other physical quantities 
show quite a rich behavior. 

Let us consider the process of an electron jumping be- 
tween sites while creating or annihilating phonons in its 
neighborhood. For example, for a chain of three sites, 
n = 3, the hopping of the electron from the first site 
will connect the states |1; 000), |2;010), |2; 100), |3; 011), 
|3; 101), |3; 020), and |3; 110). Here a vector is expressed 
as \j; ...mo, mi, ...), where j is the electron position, and 
TOfc refers to the number of phonons on site k. It is in- 
teresting to see that the connectivity of this portion of 
Hilbert spacephas the structure of a Cayley tree, or a 
Bethe lattice.Ej Under the condition of resonance, i.e., 
A = ui, all the states above for n = 3 have the same en- 
ergy. The off-diagonal matrix elements therefore break 
the degeneracy and allow level mixing, and one can ex- 
pect that some kind of band may form in the limit of large 
n. As we will see later on, and in contrast to the resonant 
case, the spectrum in off-resonance is still|-cpmposed of 
ladders although with a complex structureO 

In the basis we list above, the Hamiltonian for n = 3 
takes the form 
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A direct consequence of the Cayley tree connectivity of 
the near degenerate basis is that the Hamiltonian (Q) is 
constructed by diagonal block matrices for ever larger n 
(of sizes 1, 2, 4, 2" _1 ), and mixed by off-diagonal ma- 
trix elements proportional to A. One can see that in the 
resonant case, all diagonal elements are given by —A, 
corresponding to the degenerate manifold. Notice also 
that the off-diagonal elements are given by ±A, except 
for a few sporadic elements which appear as V2A, as- 
sociated with higher number of phonons in a site, such 
as the state |3; 020) in this case. This nearly self-similar 
structure appears for all values of n. We expect then 
that there would be little change in the physical prop- 
erties in the large n limit if we replace —y/2X in those 
few off-diagonal spots with —A. This will be confirmed 
by our numerical calculation, as we will see below. After 
the substitution, the new symmetrized Hamiltonian H sym 
takes a full self-similar form, and allows one to exactly 
solve analytically the eigenvalue problem, and better un- 
derstand the physics of the system. Much of the behavior 
of H sym remains in the actual system described by H. 



ILL SMALL HOPPING AND WEAK COUPLING 
REGIME 

Let us consider the eigenvalue problem for the Hamil- 
tonian H sym , defined as that in (^), except that the el- 
ements — v2A are replaced by —A. By using the block 
decomposition formula 

det ( C D ) = det ( D ) dct ( A - BD ~ lc )i (8) 

we find that the eigenvalues of the Hamiltonian H sym are 
determined by the equation 

e^e^ s ...si_ 1 = 0, (9) 

where N = 2", e = £ — t n , s is the energy eigenvalue, t k 
is defined by 

i fc = (n-fc)(w-A)-w, (10) 

and 

£k+1 =e-t k -^. (11) 



Then, e k can be written as a continuous fraction in k 
steps 

_ P f 2A 2 

£k+i — £ - tk 2p • 

e-tk-i ^2 (12) 

£ - tk-2 



A. Resonant regime (A = uS) 

In this case, t k = (after a constant energy shift of to 
is made), we obtain the eigenvalues of H sym as 

e fcj =2V2Acos(^ T ) , (13) 

where k = 1,2, ...,n, and j = l,...,fc. The asymptotic 
bandwidth of the spectrum is therefore given by 4v/2A. 
Notice that this behavior is similar to that of a system 
in magnetic fields. In that case, the electron-phonon in- 
teraction breaks the degeneracy (under the condition of 
magnetic resonance) and leads to a "resonance splitting." 
In that case, the splitting is oc a 2 / 3 for a three dimen- 
sional system, and oc a 1 / 2 for a two dimensional sys- 
tem, whftee a is the Frohlich electron-phonon coupling 
const ant .til 

Let us consider the question of eigenvalue degeneracy 
by solving s^ e^ e^ e^ 6 — 0, which will provide us 
with a good estimate in the large N limit. After some 
algebra, one obtains for the eigenvalue equation, 

£ 5JV / 16 (e 2 -2A 2 ) JV / 8 ( £ 2 -4A 2 ) JV / 16 
x (e 4 - 6e 2 A 2 + 4A 4 ) Ar / 16 = . (14) 

We then have eigenvalues e — 0, ±V2A, ±2A, and 
±Av 3 ± \/E, with respective degeneracies of 5 A/16, 
N/8, N/16, and N/16. Higher order polynomials would 
slightly improve the degeneracy estimate for successively 
higher eigenvalues. 

Our numerical results for the eigenvalues of H sym , Fig. 
1, and H, Fig. 2, exhibit very similar characteristics, as 
expected. The different panels in these figures show spec- 
tra for n — 8, 9, and 11. One can see that the main 
structure of the spectra changes little with increasing n, 
although finer structure appears for higher n, filling to 
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FIG. 1: Energy spectra for fully self-similar H aym for systems 
with different lattice size, (a) n = 8, (b) n = 9, and (c) n — 11. 
Notice the clear fractal structure of spectra, as anticipated 
from the structure of H syrn . Increasing n produces finer scale 
structure in the gaps. Here A = U) = 1, and A = 0.1, so that 
the asymptotic bandwidth is A\[2\ — 0.5657. 

some extent the gaps of the previous structure. Our an- 
alytical results match exactly those of Fig. 1. The self- 
similar structure of the spectrum is the manifestation of 
such symmetry in H sym . One can see that the origi- 
nal highly degenerate manifold is broadened into a semi- 
continuous band by the off-diagonal mixing elements of 
H sym . Notice however that large residual degeneracies 
remain at the center of the band, and other symmetrical 
values, as described above by Eq. ([Ti]). 

The structure of the full Hamiltonian H is shown in 
Fig. 2, which still maintains a nearly self-similar struc- 
ture, despite the sporadic — 'asymmetrical' terms. 
Notice however that the degeneracy at e = and other 
'plateau values' is not exact here, but only seen as a slight 
break (or slope) of each plateau. 

Figure |^ compares the results for the bandwidth of H 
and H S ym from both analytical and numerical results. 
Notice that 4\/2A is the asymptotic (large n) analytical 
result from Eq. (|l^), for H sym . We see in Fig. || that the 
numerical results for the bandwidth of H sym converge 
quickly to the analytical prediction. Similarly, the nu- 
merically obtained bandwidth of H is only slightly larger 
than the bandwidth of H sym , and clearly has also a fi- 
nite asymptotic value. This is an interesting feature of 
H, that although one has large Hilbert space dimension 
N - 1 = 2" - 1 (= 2047 for n=ll, for example), one 
obtains a finite bandwidth. This is of course also asso- 
ciated with the fact that there are large degeneracies in 
the energy spectra, which increase with N (see Eq. fli] ). 

To get better understanding of the physics, we show 
characteristic electron probabilities, as well as phonon 
spatial distributions, for different eigenstates in Fig. 
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FIG. 2: Energy spectra for H of systems with different lattice 
size, (a) n — 8, (b) n = 9, and (c) n = 11. Notice the 
approximate fractal structure of the spectra. Increasing n 
gives similar structure. Here A = u) = 1, A = 0.1, and the 
bandwidth slightly increases from one panel to the next. A 
and B labels in panel (a) indicate states described in Fig. [l] 
and |. 
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FIG. 3: Bandwidth results for A = 0.1 and different values of 
n. The asymptotic bandwidth for H aym is 4\/2A = 0.5657, 
shown as dashed line. We see that bandwidth for H is W > 
6A. 

and |^. The electronic probability function is given by 

Pud) = W = £ {ro} iq {m} i 2 . (is) 

where the coefficients in the last expressions are obtained 
from the diagonalization of H , so that the eigenvectors 
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are given by 



\i; {m}} , 



(16) 



for each z/-eigenstate. Similarly, the corresponding 
phonon spatial distribution is given by 



N v (j) = {* v \aj aj \* v 



— Hu{m} l^{r. 



J ! 



(17) 



for each eigenstate. These two spatial distribution func- 
tions give us an idea of how the two different components 
of each state are related to one another. These functions 
are a projection of the rather subtle coherent interactions 
(or mixtures) between the electron and phonon subsys- 
tems. One can see in Fig. 4a that for the non-degenerate 
state A in Fig. 2a, the electron is extended throughout 
the n = 8 lattice. At the same time, the phonon am- 
plitude is also extended along the lattice, and one can 
picture the phonon 'cloud' as 'surrounding' the electron 
all along in Fig. 4b, effectively describing a 'tunneling 
polaron'. In contrast, the degenerate state B at e = 
in Fig. 2a, has its electron component highly localized at 
the right end of the lattice, Fig. 5a, while the phonon 
cloud is away, Fig. 5b, and nearly 'detached' from the 
electron. One can describe this as a 'stretched polaron' 
(a precursor of the polarop-dissociation predicted at high 
fields in polymer systemsEHl). 

Although by definition the electronic distribution func- 
tion P v is normalized to unity, the phonon function N v is 
not necessarily so, since the electron can create and ab- 
sorb phonons as it propagates up and down the lattice. 
This difference in phonon content for each state can be 
seen by comparing Fig. 4b, and 5b, as the total phonon 
number is clearly larger in the latter (stretched polaron) 
case. To study this feature throughout the spectrum, 
we show in Fig. 6 the total average phonon number of 
all states in the chain n = 8. One can see clearly that 
states at the center of the band (e w 0) have many more 
phonons (« 7 = n— 1) than the rest. An estimate for the 
average number of phonons for different states can be ob- 
tained by the following argument. The average phonon 
number is given by (N v ) = P u (i)rrii. For an extended 
state, P v ~ 1/n, for a chain of n sites. A state with an 
electron at site i can be obtained by the electron hopping 
i — 1 steps from the first site, while emitting rrii rs i — 1 
phonons in the process. For extended states, this yields 
(N) ~ ^2i(i — l)/ n = i n ~ l)/2, which matches well 
the numerical results in Fig. 6, where the most extended 
states have (N) ~ 3.5 (such as the state labelled A). In 
contrast, an electron localized at site n, as B in Fig. 6, 
has emitted n — 1 phonons in the process, just as one 
finds in the exact calculation. 



1. Transport through structure 

We should emphasize that the properties illustrated 
in Fig. 4 and 5 are quite generic: states in one of the 
highly (nearly) degenerate levels show different degrees 





FIG. 4: (a) Electron probability P(j) for each site j, and 
(b) spatial distribution of phonons N(j) for state labelled A 
in Fig. 2a. Large electron amplitude and phonon number 
throughout describe an extended, delocalized polaron, jointly 
illustrated in the composite figure at bottom. 



of electron localization and an abundance of phonons, 
in a stretched polaron configuration. In contrast, states 
with non-degenerate companions are delocalized polarons 
throughout the chain with low phonon content. Even 
though is the presence of phonons which delocalizes the 
electron (in the sense of the original WS ladder), we 
obtain that the localized electrons are accompanied by 
many phonons, reminiscent of the self- localized polarons 
when the coupling is strong. 

The rather complicated behavior of the states would 
of course be reflected in various properties of the system. 
Let us focus here on what one could measure if electrons 
where injected from one end of the chain and were col- 
lected out of the other end. This is motivated by the 
ability to carry out just such experiments under strong 
electric fields in semiconductor superlattices, as well as 
in other systems. The transport properties provide infor- 
mation about the density of states in the structure and 
the spatial charge distribution of the different states (for 
optical response, sec llq]). The contribution of various 
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FIG. 6: Total average phonon number for each eigenstate 
in Fig. 2a. Higher value plateaus are associated with the 
more localized electronic amplitudes, such as state B. Lowest 
phonon counts correspond to extended electronic states, such 
as A. 




FIG. 5: (a) Electron probability P(j) for each site j, and 
(b) spatial distribution of phonons N(j) for state labelled B 
in Fig. 2a. Electron is localized at right end of the chain, 
while phonons lie throughout, yielding a stretched polaron, 
illustrated at bottom figure. 



eigenstates to the charge transport can be described by 
the quantity 



D„=P„(l)P„(n), 



(18) 



since the tunneling amplitude is proportional to the den- 
sity of states in the leads and the wavefunction ampli- 
tudes at both ends of the structure (the site j : = 1 
and j — n)x3 In fact, the tunneling probability can be 
calculated from the S'-matrix formalism as \T\ 2 , where 
T = (e f , R\S\Ei, L) . Here R and L refer to the right and 
left leads. In the wide band limit, 



T oc T J dtidhe i( - e f t2 - eit ^{n\G R (h-ti)\l) 
= ^2r{n\ V ){ V \l)5(e f - ei ), (19) 

V 

where T describes the electron interaction with contacts, 
and G R is the retarded Green function connecting both 
ends of the structure (j = 1 and j = n). We see then 



that the quantity D v = |(n|^)(^|l)| 2 oc |T| 2 is directly in- 
volved in the tunneling probability. Figure 7a shows the 
density of states and 7b shows the quantity D v for the 
system with n = 8. We can see that all the states at the 
center of the band contribute zero to the transport am- 
plitude through the chain, and this behavior is exhibited 
by all the high peaks in the DOS. This is clearly con- 
sistent with the spatially localized charge nature of the 
highly degenerate states. On the other hand, D shows 
large values for states 'in the gaps', confirming in fact 
that the non-degenerate states in the spectrum have an 
extended nature. The variations shown in D would then 
be reflected in strong amplitude modulations within each 
of the phonon replicas in tunneling experiments,!^ when- 
ever the resonance regime is reached. 



B. Non-resonant regime (A 7^ ui) 

By studying the equivalent continuous fraction ( |T2| ) for 
the non-resonant regime, we find that the spectrum for 
this case is quite different from that under resonance. 
For the case of small coupling constant, i.e., X/u> <C 1, 
and A/ A <C 1, the unperturbed spectrum is a series of 



ladder levels £j = — {w — A)j (j = 1, . 



instead of 



the degenerate manifold in the degenerate case. Notice 
that the spacing between the ladders is given by the de- 
tuning, A' = \ui — A|, instead of the original WS lad- 
der energy. The correction introduced by the interaction 
breaks some of the detuned ladder degeneracies and gives 
some substructure to the ladder. Solving the first three 
factors, e^ 2 ef^ 4 e^ 8 = 0, to second order, A 2 , we find 
that the spectrum is a series of states with energies nA', 
nA' + 2A 2 /A', (n - 1) A', (n - 1) A' - 2A 2 /A', etc., with 
respective degeneracies N/4, N/A, N/8, N/8, etc. No- 
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FIG. 7: (a) Density of states, and (b) tunneling probability 
function D„ for n = 8 system. 
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FIG. 8: Energy spectra for system with lattice size n = 8. (a) 
For full Hamiltonian H, and (b) for H sym . Here 2A = u = 1, 
and A = 0.1. A' and B' states are shown in Fig. 9. 



tice that here the electron-phonon produces level split- 
tings ~ A 2 , unlike the resonant case where the splitting 
(bandwidth) is linear in A. The problem is also solved 
numerically and we show the resulting spectra in Fig. 8, 
for both H (Fig. 8a), and H sym (Fig. 8b). Our analytical 
results match the numerical results for H sym quite well. 
We notice again that there are only small differences be- 
tween the spectra of H and H sym . 
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FIG. 9: (a) Electron probability P(j) for each site j for states 
A' and B' indicated in Fig. ^a. (b) Corresponding phonon 
distributions N(j). 



Figure 9a, shows the electron probability for different 
states in the spectra, labelled A' and B' in Fig. 8a. Figure 
9b shows the distribution of phonons for the same states. 
It is apparent that the electron wave function is localized, 
and separated from its phonon cloud. The two states only 
differ, basically, on which site to be localized about (as 
one would expect from their position in different 'rungs' 
in the ladder in Fig. 8a). In fact, the different states in 
the A' ladder are associated with different lattice sites 
(the highest rung states at e ~ nA' being localized at 
j ~ n). We also show in Fig. 10 the average phonon 
number per state for the entire spectrum. The structure 
of this figure is quite different from that in Fig. 6, as 
they reflect dramatically different dynamical behavior. A 
phonon counting argument as the electron hops and emits 
phonons can explain the average number of phonons in 
the different states in Fig. 10, although clearly here it 
produces a nearly negligible mixing of states (~ A 2 ). 

One can further explore the extent of the localization 
nature of the eigenstates. In Fig. 11, we show the elec- 
tronic probabilities of state (A') for various A' (or elec- 
tric field) values. We might expect that the localization 
length would be 2A/A', like the localization length for 
the original WS state is 2t/A. This simple expectation 
is confirmed by numerical calculations, as seen in Fig. 
12, where the solid line shows the 2A/A' dependence, 
and the dots indicate the localization length for states 
of various A' [the localization length L is extracted from 
a fit to an exponential amplitude drop about its central 
site, e - ! 2; !/ i ]. One would in fact expect that the mani- 
festation of localization in real space (and energy) would 
be susceptible to be measured via Bloch oscillations, just 
as in a typical WS ladder. In this case, however, the 
frequency of Bloch oscillations is A', and not A. Thus 



8 




100 150 200 250 300 

State index 

FIG. 10: Total phonon number for each eigenstate in the non- 
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FIG. 11: Electron probability P(j) for state A' in Fig. ^a for 
various detuning values A'. 



case the dominant conduction behavior, if any, would be 
t' " ' ' d 
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FIG. 12: Relation between localization length L and 1/A'. 
Solid line shows L/d — 2A/A' dependence, while the dots 
show value of L obtained from fit of functions in Fig. [y] to 
the form exp( — \x\/L). 



IV. BEYOND SMALL HOPPING AND WEAK 
COUPLING REGIME 

All the discussions so far have been focused on the 
Hamiltonian H (or H S y m ), valid only under the condi- 
tion of small hopping constant t/A <C 1. In this section 
we extend our discussion beyond this regime, which in 
principle may be relevant, depending on different system 
parameters. 

We use the Lang-Firsov canonical transformation 
H* = e s Hoe~ s , where S = —g^2j CjCj(aj — Oj+i), and 
g = j/uj. The transformed Hamiltonian takes the form 



the electron-phonon interaction changes the frequency of 
Bloch oscillations, and this effect may be possibly ob- 
served in experiment. 

From the discussion, it is clear that the eigenstates 
are localized in the non-resonant regime. As such, they 
contribute little to transport across the chain. One can 
easily calculate the quantity D v as before. As expected, 
the value of D u is very small (< 10~ 12 ) for every state, 
and very different from the resonant case. This behav- 
ior (not shown), is completely consistent with the local- 
ized nature of the eigenstates. Thus, we see that in this 



3 

t^ 92 X[ C J + +l C J e9(aj++1 " aj+)e_S(Qj + 1_aj) + h - C -} ■ 
3 



(20) 



One can see that this Hamiltonian is suitable for study- 
ing strong coupling dynamics. As in the case of H, we 
use the operators for WS ladder states, dj. Then the 
Hamiltonian can be rewritten in the form 



J 



5 t'J 2 (2i'/A)V-[(< 



'3 + 1 



*J+1 



■)(dtd 



3+1 



d i+i d j)} 



(21) 
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where t' = te~ 9 . We see then that the effective hopping 
constant t' becomes smaller than t. Thus the problem can 
be mapped to that in the previous section, with effective 
hopping constant t' , and effective coupling constant A', 

t' = te-9 2 , X' = gt'Jl (2t'/A) . (22) 

It is clear that A' — * 0, when g — ► 0, as expected, and 
the spectrum is the discrete WS ladder in the absence of 
electron-phonon interaction. Moreover, one can see that 
for weak coupling, the bandwidth W ~ 6A' is propor- 
tional to g, just as in the situation discussed in III. A (see 
Fig. 3). Beyond the weak coupling and small hopping 
regime, however, W and A' have a nonlinear dependence 
on g. Figure 13 shows the coupling constant dependence 
of the bandwidth W for the resonant and stronger cou- 
pling limit, A = to = t = 1. Notice that the bandwidth 
W ~ 1 for 7 > 0.6. At this coupling, the bands in 
neighboring WS rungs will begin to overlap (and inter- 
band terms would need to be considered) . It is clear that 
these corrections to the dependence of W appear also as 
function of t, A, and/or oj. 

Similar discussions are also relevant in the off-resonant 
case. For example, since the localization length is pro- 
portional to 2A'/A' (section III.B), there is also a cor- 
rection to the simple relation L oc 2A/A', when t and/or 
g are not small. For example, with increasing coupling 
constant, the correction to the ladder ~ 2A' 2 /A' may be- 
come comparable with the spacing between the detuned 
ladder levels A'. In that case, the ladder structure for the 
off-resonance case will also disappear. In this situation, 
the electron can also become delocalized, even if in the 
non-resonant regime, but due to the strong coupling 7. 

It is interesting to, compare our results to those of 
Bonca and Trugmanjj as one goes from the weak to the 
strong tunneling regime. These authors found by nu- 
merical calculation of the drift velocity that for small 
electron-phonon coupling, there are energies where the 
electron cannot propagate. In our approach, this corre- 
sponds to the formation of the quasi-continuous band for 
each of the WS ladders. The electron cannot propagate 
when its energy lies in the band gap. However, with in- 
creasing coupling constant, the band gap will disappear 
because of the overlap of subsequent bands. More quan- 
titatively, let us consider the case t = u> = 1 (as the 
parameters in Q). Since the hopping constant is large, 
one needs to adopt the formalism in this section. As seen 
in Fig. 13, the bandwidth increases with increasing cou- 
pling constant, producing band overlap for large enough 
7. This vanishing of the gaps as a function of 7 would co- 
incide with the resumption of particle drift in the system. 
Just such behavior was found numerically in , and one 
can now explain it in terms of the level structure of the 
system. It is also interesting to notice that the overall 
envelope seen in D v is similar to that of the drift veloc- 
ity: small at the center and edge of each band, as shown 
here in Fig. 7b, and qualitatively like the graphs in figure 
4 of reference . 
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FIG. 13: Bandwidth in the resonant and stronger coupling 
regime as a function of coupling constant 7. Here, A = u = 
t = 1. As 7 > 0.6, neighboring bands will overlap, indicated 
here with dotted line. 

V. CONCLUSION 

We have studied the coherent dynamics of Holstein 
polarons in a strong electric field. We have found that 
with the help of phonons, a sort of quasi-continuous band 
will form under resonance conditions, even for weak cou- 
pling constants. The band shows approximate fractal 
self-similar structure, which is inherent in the near self- 
similarity of the Hamiltonian. Although the phonons can 
help the electron jump from one WS ladder state to an- 
other, the phonons can also prevent the electron from 
propagating, if too many phonons are involved. These 
peculiar interaction differences give rise to a variety of 
unusual states, including: (a) delocalized polarons de- 
spite the strong electric field, with a typical phonon cloud 
accompanying the electron; and (b) states with high de- 
generacy at the band center, where the electron is lo- 
calized in a site of the lattice, and the phonon is located 
away from the electron, in a stretched configuration. The 
band structure is also manifested in transport properties 
of the system, which we expect could be observed in tun- 
neling experiments. The level structure and extension of 
the different states will appear as a modulation of the 
phonon replicas in the tunneling experiment. For weak 
coupling, this would only occur if the system is in res- 
onance, uj — A, as away from that condition the states 
are basically localized and would not transport current 
(except for thermal effects). In a given structure, the res- 
onance condition can be reached by sweeping the electric 
field, while monitoring the tunneling through the struc- 
ture. As the spacing between the deformed rungs and the 
localization length of the eigenfunction can be adjusted 
by electric field, it would be quite interesting to see the 
transition in experiments. With increasing coupling con- 
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stant, the 'minibands' will overlap, and give rise to an 
overall merging of the phonon replicas in tunneling, even 
when away from the resonance regime. 

Acknowledgments 

We thank A. Weichselbaum for help with graphics. 
This work was supported in part by US DOE grant no. 



DE-FG02-91ER45334, and the Condensed Matter and 
Surface Sciences Program at Ohio University. SEU ap- 
preciates the kind hospitality of K. Ensslin's group during 
his stay at ETH. 



1 E. E. Mendez, F. Agullo-Rueda, and J. M. Hong, Phys. 
Rev. Lett. 60, 2426 (1988); E. E. Mendez and G. Bastard, 
Phys. Today 46, 34 (1993). C. Waschke, H. G. Roskos, 
R. Schwedler, K. Leo, and H. Kurz, Phys. Rev. Lett. 70, 
3319 (1993). C. Martijin de Sterke, J. N. Bright, P. A. 
Krug, and T. E. Hammon, Phys. Rev. E 57, 2365 (1998); 
S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Q. Niu, 
and M. G. Raizen, Phys. Rev. Lett. 76, 4512 (1996); M. 
B. Dahan, E. Peik, J. Reichel, Y. Castin, and C. Salomon, 
Phys. Rev. Lett. 76, 4508 (1996). 

2 R. Tsu and G. H. Dohler, Phys. Rev. B 12, 680 (1975). 

3 D. H. Dunlap and V. M. Kenkre, Phys. Rev. B. 34, 3625 
(1986); Phys. Lett. A 127, 438 (1998). X.-G. Zhao, Phys. 
Lett. A 155, 299 (1991); 167, 291 (1992). M. Holthaus, 
Phys. Rev. Lett. 69, 351 (1992). W. Zhang and X.-G. Zhao, 
Physica E 9 , 667 (2001). 

4 X.-G. Zhao, R. Jahnke, and Q. Niu, Phys. Lett. A 202, 
297 (1995). 

5 A. W. Ghosh, L. Jonsson, and J. W. Wilkins, Phys. Rev. 
Lett. 85, 1084 (2000). T. Dekorsy, A. Bartels, H. Kurz, K. 
Kohler, R. Hey, and K. Ploog, Phys. Rev. Lett. 85, 1080 
(2000). 

6 D. Emin and C. F. Hart, Phys. Rev. B 36, 2530 (1987). 

7 J. Bonca and S. A. Trugman, Phys. Rev. Lett. 79, 4874 
(1997). 

8 V. V. Bryxin and Y. A. Firsov, Solid State Commun. 10, 
471 (1972). V. L. Gurevich, V. B. Pevzner, and G. Iafrate, 
Phys. Rev. Lett. 75, 1352 (1995). 

9 A. O. Govorov and M. V. Entin, JETP 77, 819 (1993). 

10 A. O. Govorov, Solid State Commun. 92, 977 (1994). 

11 L. I. Korovin and S. T. Pavlov, Sov. Phys. JETP 26, 979 



(1968); S. Das Sarma and A. Madhukar, Phy. Rev. B 22, 
2823 (1980); F. M. Peeters and J. T. Devreese, Phys. Rev. 
B, 31, 3689 (1985); 36, 4442 (1987); J.-P. Cheng, B. D. 
McCombe, J. M. Shi, F. M. Peeters, and J. T. Devreese, 
Phys. Rev. B 48, 7910 (1993). 

S. V. Rakhmanova and E. M. Conwell, Appl. Phys. Lett. 
75, 1518 (1999). 

P. M. Petroff, A. Lorke and A. Imamoglou, Phys. Today 
54, 56 (2001). 

T. Holstein, Ann. Phys. (N.Y.) 8 , 343 (1959). 

W. Zhang, A. O. Govorov, and S. E. Ulloa, Europhys. 

Lett., to appear (2002). 

For references on Cayley tre e properties, see fo r example 
K. Kundu and B. C. Gupta, |cond-mat/9705150| . 
The basis described here ignores the possibility of the 
electron hopping against the field, while creating phonons 
[absorption is accounted for in this basis when jumping 
against the field, however). Basis states such as |2; 012) or 
|2; 021), have a larger energy, produce only a small correc- 
tion of the resulting states, and are ignored here. These 
additional basis vectors can be straightforwardly included, 
although greatly increasing the relevant space and asso- 
ciated matrices. Their inclusion is more important in the 
case of strong tunneling and/or coupling constant. 
N. Zou, K. A. Chao, and Y. M. Galperin, Phys. Rev. Lett. 
71, 1756 (1993); C. A. Stafford and N. S. Wingreen, Phys. 
Rev. lett. 76, 1916 (1996). 

J. G. Chen, C. H. Yang, M. J. Yang, and R. A. Wilson, 
Phys. Rev. B 43, 4531 (1991). 



